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LONG-TERM  GOALS 

Our  long-tenn  objective  is  to  establish  a  database  of  sediment  and  bottom  properties  including 
sediment  compressional  wave  speed,  shear  wave  speed,  attenuation  and  density.  Using  this  we  propose 
to  provide  the  best  sediment  model  for  any  region  with  estimates  of  model  uncertainty.  These  data 
would  be  appropriate  for  incorporation  in  state-of-the-art  propagation  models  for  acoustic  effects 
prediction  on  the  marine  environment. 

OBJECTIVES 

1 .  Develop  a  surficial  sediment  model  for  the  world’s  oceans  based  on  optimal  estimation 
techniques,  geological  models  and  data. 

2.  Integrate  the  sediment  model  into  WHOI  oceanographic  models  for  input  to  ray/  normal  mode/ 
PE  based  sonar  effects  studies. 

3.  Provide  estimates  for  model  uncertainty. 

4.  Apply  this  model  in  selected  regions  to  test  its  accuracy  and  effectiveness. 

APPROACH 

As  a  first  step  we  plan  to  collect  all  the  available  sediment  property  data  at  the  region  of  interest.  This 
includes  all  the  geoacoustic  (compressional  wave  speed,  shear  wave  speed,  attenuation  etc.)  and 
geotechnical  properties  (bulk  density,  permeability,  porosity  etc.).  The  geotechnical  properties  can  be 
used  to  derive  geoacoustic  properties  using  a  suitable  model.  These  data  can  be  obtained  from 
literature  or  from  available  sediment  property  databases.  In  addition  to  these  direct  measurements  of 
sediment  properties  we  intend  to  collect  other  parameters  which  influence  sediment  distribution  at  any 
given  region.  These  factors  include,  water  depth  at  the  location,  distance  to  the  shore,  region,  latitude, 
current  speed  and  proximity  to  major  river  discharge 

After  collecting  the  data  as  described  above,  the  next  step  is  to  express  the  non-linear  model-data 
relationship  in  the  form: 

d=G  (m) 

where 


1 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
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The  matrix  G  defines  the  non-linear  relationship  between  the  sediment  data  (m)  and  the  factors  which 
influence  them  (d).  The  main  task  then  will  be  to  generate  the  kernel  G  based  on  the  relationship 
between  d  and  m.  Various  sediment  geoacoustic  models  (Hamilton4,  Bowles5  etc.)  sediment  transport 
models  and  geological  information  can  be  used  to  develop  the  kernel  matrix  G.  We  plan  to  proceed  as 
follows: 

1 .  Collect  all  the  geoacoustic  parameters  available  at  various  sites  in  the  region  of  interest. 

2.  Collect  also  other  parameters  of  influence  like  water  depth,  latitude,  distance  from  shore,  etc. 

3.  Develop  a  relationship  (G)  between  these  parameters  using  half  of  the  available  data. 

4.  Using  G,  predict  the  remaining  half  of  the  data  and  compare  with  the  true  data  and  thereby  test 
the  model.  Tune  the  model  until  there  is  good  agreement. 

5.  Using  G  calculate  the  geoacoustic  parameters  at  any  location. 

WORK  COMPLETED 

A  sediment  model  suitable  for  the  Southern  California  Bight  (SOCAL)  region  was  first  developed.  The 
second  test  case  corresponds  to  the  Middle  Atlantic  Bight  region.  A  program  was  written  to  generate 
sediment  parameters  for  the  Middle  Atlantic  Bight  region  which  was  chosen  as  the  test  site  in  the  last 
(August  2001)  ESME  workshop.  The  location  lies  between  latitudes  37  N  and  43  N  and  66  W  and  74 
W.  The  program  is  written  in  such  a  way  that  given  the  location  and  the  frequency  of  interest,  the 
sediment  properties  are  estimated  using  all  the  available  data.  The  sediment  information  available 
includes  core  data,  sediment  maps,  sediment  type  information  etc.  For  this  location  we  were  able  to 
find  some  historical  data  which  is  listed  below. 
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1 .  Core  data: 

19 

♦  Atlantic  Margin  Coring  Sites  6009  to  6021  ’  .  The  sediment  type  information,  sediment 
properties,  layer  thickness  etc.  were  obtained  from  these. 

♦  Gravity  core  data  .  Gravity  cores  were  taken  and  analyzed  from  five  locations  as  part  of 
Primer  study. 

♦  Vibrocore  data 

♦  Piston  Core  data  from  USGS  mapping  program4 

2.  Sediment  type  information:  from  USGS  sediment  maps  for  the  area5. 

3.  Compressional  wave  speeds  and  layer  thickness  using  high  resolution  seismic  studies  (Brocher  and 
Ewing) 

4.  Sediment  sound  velocity  functions  for  North  Atlantic  (Houtz) 

5.  Sediment  tomography  using  Primer  data  (Potty  et.al.) 

6.  Thickness  of  surficial  sand  sheet  (Knebel  and  Spiker) 

7.  Compressional  wave  velocities  from  multichannel  refraction  arrivals  (Me  Ginnis  and  Otis) 

8.  Compressional  and  shear  wave  speeds  using  gravity  wave  inversion  technique  (Trevorrow  and 
Yamamoto) 

9.  Sediment  and  geoacoustic  models:  Used  for  predicting  the  compressional  and  shear  wave  speeds  as 
well  as  attenuation  as  function  of  frequency  for  different  sediment  types.  These  include  Hamilton 
model1 1 ,  Bowles12,  and  Biot-  Stoll  model13. 

o 

10.  Sediment  layer  thickness  maps  :  Used  for  predicting  the  sediment  thickness. 


Sediment 


Compr.  Speed,  Shear  Speed 
Attenuation,  Density 


Basement 


Compr.  Speed,  Shear  Speed 
Attenuation,  Density 


Figure  1.  Two  layer  model  to  represent  the  sediment  properties  used  in  this  study. 

The  program  uses  these  data  to  arrive  at  the  sediment  parameters  at  a  given  location.  The  sediment 
compressional  speed  is  estimated  using  the  core  data,  sediment  type  infonnation  and  the  Hamilton’s 
values  for  these  sediment  types.  Sediment  thickness  is  estimated  based  on  the  core  data  and  sediment 
thickness  map.  Shear  speed,  and  Compressional  and  shear  attenuation  are  determined  mainly  using  the 
sediment  type  infonnation  and  predictions  by  Hamilton  and  Stoll.  The  bottom  properties  are  specified 
for  a  simple  two-layer  model.  Some  of  the  properties  (compressional  and  shear  speeds  in  the  basement) 
are  assumed  constant  throughout  the  area.  The  output  parameters  are  written  in  a  text  file. 
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Figure  2.  Left  panel  shows  the  location  of  the  ESME  test  case.  Right  panel  shows  the  sediment 
sound  speed  obtained  using  the  sediment  model  developed. 


IMPACT/APPLICATIONS 

Our  effort  is  intended  to  provide  the  best  possible  estimates  of  sediment  data  needed  for  the 
propagation  modeling  component  of  the  ESME  initiative.  This  will  enable  the  propagation  modelers  to 
use  ‘best  guesses’  when  direct  estimates  of  the  sediment  data  is  not  readily  available. 

TRANSITIONS 

We  expect  that  this  sediment  model,  when  complete,  will  be  useful  to  the  acoustic  community  as  a 
whole  as  an  important  database  of  sediment  property  information.  This  model  can  be  integrated  into 
OAML  based  on  discussions  with  NAVOCEANO  personnel. 
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